REVIEW ARTICLE 



Spin-orbit coupling and electron spin resonance for 
interacting electrons in carbon nanotubes 

A De Martinof , R Eggerf §, F Murphy- Armando! an d K 
Hallberg| 

f Institut fur Theoretische Physik, Heinrich-Heine-Universitat, D-40225 Diisseldorf, 
Germany 

| Institute Balseiro, Centra Atomico Bariloche, Comision Nacional de Energi'a 
Atomica, 8400 S.C. de Bariloche, Argentina 

Abstract. We review the theoretical description of spin-orbit scattering and electron 
spin resonance in carbon nanotubes. Particular emphasis is laid on the effects of 
electron-electron interactions. The spin-orbit coupling is derived, and the resulting 
ESR spectrum is analyzed both using the effective low-energy field theory and 
numerical studies of finite-size Hubbard chains and two- leg Hubbard ladders. For 
single-wall tubes, the field theoretical description predicts a double peak spectrum 
linked to the existence of spin-charge separation. The numerical analysis basically 
confirms this picture, but also predicts additional features in finite-size samples. 



PACS numbers: 71.10.-w, 73.63.Fg, 76.30.-v 



§ To whom correspondence should be addressed (egger@thphy.uni-duesseldorf.de) 



Interacting electrons in nanotubes 



2 



1. Introduction 

Nanotubes constitute a new class of mesoscopic quantum wires characterized by the 
interplay of strong electron-electron interactions, low dimensionality, disorder, and 
unconventional spin dynamics [Il|2llI3llllElE10IE|- I n a sense , they represent an ideal 
model for strongly correlated mesoscopic systems, where in fact basically all known 
effects in mesoscopic physics have been experimentally observed. Two main classes of 
nanotubes can be distinguished, namely single-wall nanotubes (SWNTs) which consist 
of just one wrapped-up graphene sheet with radius R in the nanometer regime, and 
multiwall nanotubes which contain additional inner shells [5]. Here we will focus on the 
conceptually simplest case of metallic SWNTs, where interactions should completely 
destroy the Fermi liquid picture and imply a so-called Luttinger liquid (LL) state of 
matter [HllZj. The Luttinger liquid is the generic low-energy description of metallic ID 
(single-channel) quantum wires [9J. 

Evidence for the LL behaviour of interacting ID electrons has been reported for 
charge transport in SWNTs [21 El II]- However, in such materials one also expects to 
find more dramatic consequences of the breakdown of Fermi liquid theory, most notably 
the phenomenon of spin-charge separation. This many-body effect asserts that electrons 
brought into a LL effectively break up into a charge and a spin part that travel with 
different velocities and hence will be spatially separated after some time. A recent 
proposal to detect evidence for spin-charge separation in SWNTs has been based on 
spin transport |Hj. A different (and perhaps easier to realize) proposal based on electron 
spin resonance (ESR) is reviewed in this paper, expanding on our short paper [TU]. ESR 
is a valuable experimental tool to probe the intrinsic spin dynamics of many systems. 
In ESR experiments one applies a static magnetic field and measures the absorption of 
radiation polarized perpendicular to the field direction. In the absence of SU(2) spin 
symmetry breaking terms in the Hamiltonian, the absorption intensity is then simply a 
5-peak at the Zeeman energy [TU I12j. 

Since spin-orbit (SO) interactions are generally the leading terms breaking the 
SU(2) invariance, deviations in the ESR intensity from the 5-peak, e.g. shifts or 
broadenings, are directly connected to these couplings. Below we theoretically address 
the spin-orbit interaction and the resulting ESR spectrum for interacting SWNTs, using 
both a continuum field theory and a Hubbard model description. Within the effective 
field theory, the single 5-peak is split into two narrow peaks in SWNTs if spin-charge 
separation is realized. Otherwise the ESR spectrum would form a broad band with 
thresholds at the lower and upper edge [13] • This qualitative difference is caused by 
the fact that the SO interaction in SWNTs does not spoil spin-charge separation to 
leading order. Experimental observation of the peak splitting could therefore provide 
evidence for the elusive phenomenon of spin-charge separation j^j. To experimentally 
check the predictions made below, samples free of magnetic impurities have to be used. 
Such impurities have probably spoiled previous ESR measurements for nanotubes [5 . 

The outline of this article is as follows. In Section |21 the Luttinger liquid theory of 
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carbon nanotubes is reviewed in detail, and we give an introduction to ESR theory as 
relevant for our purposes. The spin-orbit interaction is derived in Section El followed 
by a detailed discussion of the low-energy theory predictions for the ESR spectrum in 
Section HJ An alternative approach is to use numerical methods to compute the ESR 
spectrum for microscopic lattice fermion models. We shall use both a Hubbard chain 
and the more realistic two-leg Hubbard ladder formulation of interacting SWNTs [14J. 
While charge transport does not allow for such a description due to the importance of 
long-range interactions, it turns out that for ESR spectra, only short-range interactions 
are important. These are correctly captured by Hubbard-type models, and therefore 
such models are expected to be appropriate for the quantitative description of ESR 
spectra in SWNTs. We discuss this approach in Section El and compare the numerical 
results with the field-theoretical predictions. Finally, in Section |H1 conclusions and a 
brief outlook are provided. In most of the paper, we use the conventions % = c = 1 to 
simplify notation. 

2. Basics 

2.1. Luttinger liquid theory 

Starting from a microscopic lattice description of the SWNT, inclusion of the interactions 
among electrons leads to rather complicated models. In the case of short-ranged 
interactions, one can study Hubbard-type models, and we will do so in Section El For 
low-energy phenomena such as electron spin resonance, however, only bands close to the 
Fermi energy do matter. The relevant electronic properties of (not too thin) SWNTs 
are then caused only by the ir electrons of the wrapped graphene sheet. On energy 
scales \E\ < D = Tivf/R ~ 1 eV around the Fermi energy Ep (here Vp is the Fermi 
velocity), the graphene bandstructure takes a simple form allowing to develop a powerful 
field-theoretic framework for SWNTs j§| reviewed below. We mention in passing that it 
can be explicitly demonstrated that bands sufficiently far away from the Fermi surface 
will not change the results obtained from the field theory; for an explicit discussion, see 
Ref. 0. 

Simple tight-binding bandstructure calculations for graphene reveal that there are 
only two linearly independent Fermi points ("flavours") with coordinates aK in the 
first Brillouin zone (a = ±), instead of a continuous Fermi surface [I]. For \E\ < D, 
the dispersion relation around the Fermi points is highly linear (two-dimensional light 
cone). Since the basis of the graphene honeycomb lattice contains two atoms, there are 
two sublattices p — ±, and hence two degenerate Bloch states 

<P P a(r) = (2ttR)- 1/2 exp(-iaKr) (1) 

at each Fermi point a = ±. Here f — (x, y) lives on the sublattice p under consideration, 
and the correct normalization for nanotubes has been anticipated. The Bloch functions 
are defined separately on each sublattice such that they vanish on the other. One can 
then expand the electron operator in terms of these Bloch functions. The resulting 
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effective low-energy theory of graphene is the 2D massless Dirac Hamiltonian, as follows 
also from standard k ■ p theory. 

Wrapping the graphene sheet onto a cylinder then leads to transverse momentum 
quantization, and hence to the effectively ID bandstructure of a metallic SWNT. Taking 
the x-axis along the tube direction and the circumferential variable as < y < 2ttR, 
quantization of transverse motion now allows for a contribution oc exp(imy / R) to the 
wavefunction. However, excitation of angular momentum states other than m = costs 
a huge energy of the order D. Assuming that the SWNT is not excessively doped, 
following our above remarks, in the field theory we may then omit all transport bands 
except m — 0. The theory will then apply on energy scales \E\ < D and lengthscales 
larger than the graphene lattice spacing a « 0.246 nm. 

Evidently, a SWNT then forms a truly ID quantum wire with only two spin- 
degenerate bands intersecting the Fermi energy. To take this into account, the electron 
operator for spin o = ± is written as 

V a (x, V) =J2 ( Ppa{.X, y) 1ppaa(x), (2) 
pa 

which introduces ID fermion operators ip paa {x). Neglecting interactions for the moment, 
the Hamiltonian is 

#0 = -VfY,P / dx ^paadx^-paai (3) 
pacr 

which is equivalent to a massless ID Dirac Hamiltonian, where Vp = 8 x 10 5 m/sec. 

Next we discuss interactions mediated by the (possibly externally screened) 
Coulomb potential U(f — r*). The precise form of this potential will of course depend 
on details of the setup. In the simplest case, bound electrons and the effects of an 
insulating substrate are described by a dielectric constant k, and for the long-range 
Coulomb interaction, 

U(f- f ) = e2,K (4) 

y](x - x') 2 + AR 2 sin 2 [(?/ - y')/2R] + a 2 

where a z ~ a denotes the average distance between a 2p z electron and the nucleus, 
i.e. the "thickness" of the graphene sheet. Electron-electron interactions are then 
described by 




which is brought into a ID form by inserting equation (0) for the electron field operator, 
allowing to employ the large arsenal of theoretical methods available only in ID jU]. The 
result is 

H i = \Y, E J dxdx> V^ } (x - x>) (6) 

pp'cra' { ai } 

X ^la 1 cr( X )'Pl'a 2 a'( X ')%'a 3 a'( X ')%a 4 a( X )^ 
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with ID interaction potentials 

V& } (x-x') = J dydyV pai (^;, a2 (r<)U(r-f<+p^ 

These potentials only depend on x — x' and on the ID fermion quantum numbers. 
For interactions involving different sublattices p ^ p' for r and r 1 in equation (jHJ), a 
sublattice shift vector d arises 7\. To simplify the resulting ID interaction (jUJ), one 
can exploit momentum conservation. Provided we stay away from the charge neutrality 
point Ep = 0, Umklapp electron-electron scattering processes can be ignored, and the 
situation simplifies considerably. We then have only "forward scattering" processes [Zj, 
where «i = a 4 and a 2 = a 3 , and "backscattering" processes with a\ = — a 2 = a 3 = — a 4 . 
Next we introduce the potential 

r/r >\ r 2wR d y f 2nR d y' ^ ^ 

V Q [x — x = / / U r — r . 8 

For the unscreened Coulomb interaction (0J), this can be explicitly evaluated [Zj. 
From equations (JTj) and (JTJ), the forward scattering interaction potential reads V (x) + 
8p-p'8Vp(%), with 

6V p (x) = j^ R ^^P( x + P d *, V-y' + P d y) - V - V% ( 9 ) 

which is only present if r and ? are located on different sublattices. Thereby information 
about the discrete nature of the graphene network has been kept despite the low-energy 
continuum approximation. Since V (x) treats both sublattices on equal footing, the 
resulting forward scattering interaction part couples only the total ID electron densities, 

Hf ] = -j dxdx' p{x)V {x - x')p(x'), (10) 

where p(x) = Y^paa^laa^paa- This part of the electron-electron interaction is the most 
important one and is responsible for the LL behavior. Note that it is due to the long- 
ranged tail of the Coulomb interaction. The remaining interactions originate from 
short-ranged interaction processes, and since these are effectively averaged over the 
tube circumference, their amplitude is quite small, scaling as 1/R. Such couplings are 
seen below to cause only exponentially small gaps. 

For |x| ^> a, the term SV p (x) is extremely small. However, for \x\ < a, an additional 
term beyond equation arises due to the hard core of the Coulomb interaction. At 
such small length scales, the difference between inter- and intra-sublattice interactions 
matters, and SV p (0) must be computed microscopically, leading to a coupling constant 
/ characterizing the additional forward scattering contribution 

= ~f J dx iiaA-pcc'A-po.'a'^paa: (H) 

paa'aa' 

where f/a = / yfe 2 /R with a dimensionless constant 7/ depending on the tube chirality. 
An estimate for armchair SWNTs yields 7/ ~ 0.05, implying that / is very small. A 
similar reasoning applies to the backscattering contributions ot\ = —cti = 0:3 = —a^ in 
equation EJ). Because of a rapidly oscillating phase factor, again the only contribution 



Interacting electrons in nanotubes 



6 



comes from \x — x'\ < a. Furthermore, only the part of the interaction which does not 
distinguish among the sublattices is relevant, leading to 

H P = bjdx ^ia^i>-aa>%>aa>%-a*- ( 12 ) 

pp'acra' 

For the unscreened interaction (J3J, b/a = ^e 2 /R with 7^ « 7/. For externally screened 
Coulomb interaction, however, one may have b^> f . 

Further progress can be made by using the Abelian bosonization approach [Oj. 
For that purpose, one brings the non-interacting Hamiltonian (jSJ) into the conventional 
form of the ID Dirac model by switching to right- and left-movers (r = ±) which are 
linear combinations of the sublattice states p — ±. In this representation, an Abelian 
bosonization formula [3 El applies with four bosonic phase fields 6 a (x) and their 
canonical momenta U a (x). The four channels are obtained from combining charge and 
spin degrees of freedom as well as symmetric and antisymmetric linear combinations of 
the two Fermi points, a = c+, c— , s+, s— . The bosonized expressions for H and Hf ] 
read 

^0 = Ey/^[ n a + ^ 2 (^a) 2 ] (13) 

Hf ] = - f dxdx' d x e c+ (x)V (x - x')d xl 6 c+ (x'). (14) 

7T J 

The bosonized form of Hj 1 ' 2 ^ jH] leads to nonlinearities in the 9 a fields for a 7^ c+. 
Although bosonization of equation (jSJ) gives K a = 1 in equation (|TH|) . interactions will 
renormalize these parameters. In particular, in the long-wavelength limit, Hj ^ can be 
incorporated into Ho by putting 

K c+ = {l + 4:V (k ~ 0)/7rhv F y 1/2 < 1, (15) 

while for all other channels, the coupling constant / gives rise to a tiny renormalization, 
K aj L c+ = 1 + f/irhvF ~ 1. The plasmon velocities of the four modes are v a = Vp/K a , and 
hence the charged (c+) mode propagates with higher velocity than the three neutral 
modes. The dimensionless Luttinger parameters measure the correlation strength in 
the system, with the noninteracting point at K a = 1 and repulsive interactions leading 
to K c+ < 1. For the long-ranged interaction the logarithmic singularity in Vo(k) 
requires the infrared cutoff k = 2tt/L due to the finite length L of the SWNT, resulting 
in: 

K c+ = { 1 + \n(L/2nR) ) . (16) 

[ TTKhVp J 

Since Tic/e 2 ~ 137, we estimate e 2 /TiVF = (e 2 /Tic)(c/vf) ~ 2.7, and therefore K c+ is 
typically in the range 0.2 to 0.3. This estimate does only logarithmically depend on L 
and R, and therefore is expected to be almost independent of the sample under study. 
The Luttinger parameter ()16|) can also be written in the form 

2E C ^ 



K c+ = [ 1 + A 
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where E c is the charging energy and A the single-particle level spacing. The small value 
predicted here implies that a metallic SWNT should be a strongly correlated system 
displaying pronounced non-Fermi liquid effects. 

It is clear from equations (|T5j) and (fTlj) that for / = b = 0, a SWNT constitutes 
a realization of the LL. We therefore have to address the effect of the nonlinear 
terms associated with the coupling constants / and b. This can be done by means 
of the renormalization group approach. Together with a solution via Majorana 
refermionization, this procedure allows for the complete characterization of the non- 
Fermi-liquid ground state of a clean nanotube [HJ E] • From this analysis, one finds that 
for temperatures above the exponentially small energy gap 



induced by backscattering processes, the SWNT is adequately described by the LL 
model, and H) ' can effectively be neglected. A rough order-of-magnitude estimate 
is TJ, pa 0.1 niK. In the remainder, we focus on temperatures well above Tj,, where the 
nonlinearities can be neglected and the Luttinger picture applies. 

2.2. Sugawara formulation 

In the following discussion of ESR theory, it is mandatory to keep the ££7(2) spin 
symmetry explicit at all stages. To do so it is advantageous to avoid the Abelian 
bosonization used above, which breaks the spin symmetry by hand, but rather employ 
the Sugawara formulation |Hj which manifestly respects SU(2) spin invariance. This 
formulation is in fact fully equivalent to a Wess-Zumino-Witten theory, even with the 
flavour degeneracy due to the two Fermi points 0. Since our main interest is on spin 
properties, we shall however suppress the flavour index for most of what follows, but 
return to the complexities added by it later, see Section Instead of four channels, 
we then have a two-channel Luttinger liquid, with interaction constants for charge, 
K c = (K c _ + K c+ )/2 ~ (1 + K c+ )/2, and spin, K s = (K s _ + K s+ )/2 ~ 1. In Abelian 
bosonization, in order to impose the correct spin symmetry, one then fixes K s = 1 [Oj. 
This procedure is not necessary in the Sugawara treatment below. Inclusion of magnetic 
Zeeman fields B (orbital effects play no role here) only affects the spin sector and in 
general could renormalize K s . However, this renormalization is irrelevant to ESR, which 
probes the finite energy scale ~ B [TT|. In addition, the velocities v c / s = Vp/K c / S are 
used in the following. 

The Sugawara formulation uses the charge (Jl/r) an d spin currents Jr/l for right- 
and left- (r = ± = R/L) moving electrons described by chiral ID fermion operators 
ip r = R / L (x,t), where we suppress the additional spin index a. In terms of these field 
operators, chiral £77(2) spin current operators are given by 



where the colons denote normal-ordering and Pauli matrices a act in spin space. They 



ksTb = D exp[— irfivp/ V2b] 



(17) 




(18) 
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obey Kac- Moody commutation relations (//, v — x, y, z) 

[•#/*(*). J U l/r{x')\ = ±i*'(s - x')8^/A-k + ie^J x L/R {x)5{x - x'). (19) 
Likewise, charge current operators are defined as 

Jr,l =■ i>R /L (x)if> R/L (x) :, (20) 

where spin indices are summed over. Using these current operators, the Luttinger liquid 
in the £7(1) x ££7(2) invariant Sugawara formulation reads j^j, 

H = H C + H S , (21) 

with decoupled charge and spin parts. The charge sector is described by the £7(1) 
invariant Hamiltonian 

Hc = ~Y J dx ^ JrJr + JlJl : +9c : JlJr :) ' (22) 

The coupling g c is determined by the Luttinger parameter K c . Explicit expressions can 
be found in Ref. [Oj, but they are not required below. In particular, in the noninteracting 
limit, g c = and v c = Vp- The £77(2) invariant spin Hamiltonian H s commuting with 
H c is 

Itxv r -»-»-»-* r -» -» 

H s = -^r 1 dx : J R - J r + J l ■ J L : +g s dx: J R - J L :, (23) 



3 

where v s ~ is the spin velocity and g s = 6f s is the coupling constant for electron- 
electron backscattering processes. In a dynamical spin- sensitive ESR measurement, one 
should not simply discard this coupling despite the above thermodynamic argument 
invoking the smallness of gaps The LL Hamiltonian (}2~Tj) completely decouples when 
expressed in terms of spin and charge currents (or the spin/charge bosons of Abelian 
bosonization). This remarkable fact leads to the phenomenon of spin-charge separation, 
and unless the spin-orbit interaction couples spin and charge sectors, ESR will only 
probe the spin sector. 

2.3. Electron spin resonance 

Adopting the conventional Faraday configuration, the ESR intensity at frequency uj is 
proportional to the Fourier transform of the transverse spin-spin correlation function 

nn, 

= fdte iuJt (S + (t)S-(0)), (24) 



where the static magnetic field points along the z-axis, £ = J2i £« is the total spin 
operator, and £ = S x ± iS y . The Hamiltonian can be written as H = H + Hz + H\ 
where Hq represents the ££7(2) invariant nanotube model (|21|) including electron- 
electron interactions, H z = ~g e \iRH ■ S is the Zeeman term (below, often g e \iR = 1), 
and H' represents ££7(2) spin-symmetry breaking terms, in particular the SO coupling. 
Inserting a complete set of eigenstates \a) of H in equation pijl . the ESR intensity 



follows as 



m = ^J2e- Eb/kBT 5(u - (E a - E b )) \(a\S~\b)\ 2 . (25) 



a.b 
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In the absence of 777, there are only contributions from matrix elements between 
eigenstates with equal total spin S a = SV Then all states with S* = — 1 will 
lead to a 5-peak at frequency uo = B. For instance, at zero temperature, the application 
of a magnetic field B, taken as large enough to overcome a spin gap possibly present 
at B — 0, leads to a ground state with finite magnetization, So ^ 0, and the states 
with = Sq — 1 again yield the 5-peak. This can be made explicit as follows. Since 
[H + H Z ,S~] = BS~, one has (77 + H z )S-\0) = S~(H + H z )\0) + BS~\0) and 
identifying S~\0) = \a), one gets E a = E + B and thus I(uS) = Iq5(uj — B). Any 
perturbation preserving £77(2) invariance will neither shift nor broaden this peak, even 
at finite temperature To get nontrivial ESR spectra, one has to identify the leading 
perturbation breaking 577(2) invariance, which will cause finite linewidth and shift of 
the ESR peak. In metallic systems like SWNTs, one has to consider the spin-orbit (SO) 
interaction. This is done in the next section. 

3. Spin-Orbit Coupling 

3.1. Microscopic derivation 

In our derivation of the SO term, we shall neglect electron-electron interactions. Local 
electric fields exerted by other electrons on a given electron are typically weak compared 
to the ionic fields ^H], and will generally only weakly renormalize the SO couplings from 
their noninteracting values. In a single-particle picture, the SO interaction then appears 
because an electron moving in the electrostatic potential $(r) experiences an effective 
magnetic field v x V$ in its rest frame. With p = mv the SO interaction reads in 
second-quantized form 



This represents the starting point for our discussion of the SO coupling. 

For a microscopic lattice description, the electron field operator ^ a {r) [whose low- 
energy expansion is given above in equation (J2J)] can now be expressed in terms of the 
electron operators q for honeycomb lattice site i at f^, 



where xi? — fi) is the corresponding Wannier wavefunction centered at lattice site 
r*j. These localized Wannier orbitals can be chosen as real- valued functions even 
when hybridization with s-orbitals is important. For the simplest case, 2p z orbital 
wavefunctions could be used, with z perpendicular to the graphene plane, 



where Oq = h 2 /6me 2 is the effective Bohr radius. 

Inserting the expansion (|2Tjl into (H|, the SO interaction reads [Tnl M\ ITTj . 




(26) 




(27) 
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J dr$(r) [Vx(r — r}) x V%(f - f fe )] . (29) 



which explicitly breaks SU(2) symmetry. The SO vector ■u^ = —Ukj has real-valued 
entries and can be written as 

-. _ 9e^B 

4m 

The on-site term (j = k) is identically zero, and since the overlap decreases exponentially 
with \ fj — r^l, we keep only nearest-neighbour terms in equation (|28|) . Equation (|28|) 
also allows to incorporate electric fields due to impurities or close-by gate electrodes. 
To connect the SO vector to the experimentally measurable SO relaxation rate, t^q, we 
estimate this rate using Fermi's golden rule. Assuming a constant u, the probability for 
a transition from the initial state \i) = |&, t) to the final state |/) = \k', j) is given by 
Ti-tf = 2n\(f\H'\i)\ 2 5(Ei — Ef), and summing over all final states, 

r^ = 8(u 2 x + u 2 y )/3v s . (30) 

For SWNTs, SO couplings are generally expected to be small. This is in accordance with 
our approach since the SO vector ()29)1 vanishes by symmetry for an ideal 2D honeycomb 
lattice. A finite nearest-neighbour SO coupling can only arise due to the finite curvature, 
stray fields from nearby gates, or due to impurities, which break the high symmetry. 
Focusing for clarity on the curvature-induced SO coupling, the SO vector for achiral 
tubes can be seen to only depend on bond direction, 

where the nearest-neighbour bonds 5 n {n = 1,2,3) connect the two sublattices of the 
graphene sheet in real space. For explicit representations of the vectors S n , see e.g. |17j . 

3.2. Low- energy form 

From now on we then restrict our analysis to achiral tubes and construct the low-energy 
field theory Hamiltonian H' describing the spin-orbit interaction (|2"%|) . including the 
flavour indices for the moment. Using equation (J3"T|) . the microscopic form is 

3 

H ' = 1212 { ic l,p=+Un ■ <?C r - +5niP= _ + H.C.) , 
r n=l 

where r runs over all lattice sites of the sublattice p — +, and we have made the 
sublattice dependence of the lattice electron operators explicit (again spin indices are 
implicit). The expansion (j2J) implies 



which yields after ^/-integration and Taylor expansion of ip p= _^ acr (x + e x ■ 5 n ) the result 
H'= [dxJ2 e- iaR ^4 =+ , a (x) (iu n ■ a) (32) 



x 



+ H.c. 



1pp=-,a{x) + (e x ■ 0~n)d x 1pp=- ta {x) 

where we use the unit vector e x along the tube axis. Terms involving higher order 
derivatives of the field operators are highly irrelevant and can safely be ignored. In 
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equation ([32]). we have also neglected oscillatory terms that vanish in doped SWNTs 
because of axial momentum conservation. These oscillations are governed by the 
wavevector 2\/j\/vf corresponding to the doping level /i. In practice, intrinsic doping is 
unavoidable, but for a commensurate situation (// = 0), additional terms do arise that 
are ignored in this section. 

Next we define vectors given in terms of the vectors specified in (}3*Tj) . 



3 

u la = exp(-iaK ■ S n )u n 



n=l 
3 



u 2a = expf-iaif • 5 n )u n (e x ■ 5 n ). 

n=l 

We then obtain from equation (}3*2*|) the form H' = Hi + H 2 with 

H 1 = J dxJ2 ^l=+A x ) (^1,0 ' ^) ^p=-A x ) + H - c - ( 33 ) 



and 



H 2 = J dxJ2^l=+A x ) (^2,a • d x ij)p=- }a (x) + H.c. (34) 



It is then apparent that i?i will give the dominant contribution (relevant in the 
renormalization group sense, scaling dimension 1), while H 2 is marginal (scaling 
dimension 2). Since SO interactions are always small, the renormalization group 
arguments around the Luttinger fixed point with no SO couplings can be safely used 
here. Note that all contributions neglected here either violate momentum conservation 
or have scaling dimensions > 3, i.e. are highly irrelevant. 

The final step is to move from the sublattice description to the r = R/L basis using 
a unitary transformation. Noting that the term involving Reu lQ involves oscillatory 
contributions that can be neglected by virtue of momentum conservation, the dominant 
low-energy SO contribution takes the form 

H 1 = / dx^\ a - (J Lta - J Ri(X ), (35) 

J a 

with effective SO vectors 

X a = 2 Imu 1)Q = 2 Im ^ e~ ia ^ n u n . (36) 

n 

In the analysis of Sec. IH for clarity, we shall ignore the flavour index, but it is clear 
from the above how to add it to the theory. The subleading term H 2 reads (ignoring 
the flavour index, and again omitting oscillatory terms) 

H 2 = fdx V^A' -ffid^r + B.c. (37) 

r=R/L 

where 

n 
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4. ESR spectrum from field theory 

In the following, we will keep only the leading contribution ()35)) and neglect the marginal 
term H 2 . In the numerical analysis of Section the full SO interaction will be 
considered to assess the accuracy of these approximations. Since SO couplings are 
typically very small, however, the above perturbative reasoning, suggesting to drop H 2 
and the oscillatory terms, is expected to work. Remarkably, the SO interaction then acts 
exclusively in the spin sector and hence does not spoil spin-charge separation. Therefore 
the ESR intensity can be computed from the spin sector alone. 
For convenience, we define 

K=r/l=± — B ± A, (38) 

which represents an effective magnetic field acting separately on the chiral (right/left- 
moving) spin currents. The presence of two pseudo-fields already hints at the outcome of 
our calculation below, namely the emergence of a double peak spectrum. Their position 
will then be given by the absolute values A r = |A r |. The spin Hamiltonian governing 
the ESR spectrum I{ui) is then 

H = H s -J2j dxX r -J r , (39) 



with H s given in equation (J23J). The Heisenberg equation of motion for the chiral spin 
current operators is 

(d t ± v s /3d x )J r = T{g s /^)d x J- r + 9 s Jr x J L - \ r x J r . (40) 

Taking B = Be z , and using S = Jr + Jl, the ESR spectral density is 

I(u) = J dtdxe iuJt J2(J?(x>t)J^ (0,0)}. (41) 

rr' 

We then compute the ESR spectrum using the Sugawara spin Hamiltonian (|3*9*)l. An 
independent route would proceed via fermionization For simplicity, we discuss the 
case g s = (no backscattering) , but the method is flexible enough to allow treatment 
of the general case as well. However, the numerical analysis of Section El indicates that 
indeed the backscattering interaction does not affect the results in a significant way. 

Since the spin Hamiltonian decouples in the chiral spin currents, we just need to add 
the contributions due to the two chiralities. In equilibrium, using the spin susceptibility 
1/4-KVs of a chiral fermion, we have 

(/«) = JL-X* (42) 

Using the fluctuation-dissipation theorem, we can now express I(u) in ({IT} in terms of 
the imaginary part of the susceptibility tensor tt>), 

/M = £ .. £ .. 1-JL/k.T) Im * rfe - a (43) 



=± fJ.,u=x,y 
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The susceptibility can in turn directly be obtained from the equation of motion (J4(J|) . 
Formally defining the matrix (A r .) Atl/ = e M ^ a A", we find 

Xr{q ->0,u) = 7^— [iw + A r ] _1 Ar., 

47Tf s 

and straightforward algebra then yields the ESR spectrum 

'M = E [1 + W 4 ^ (1 _ A ;_ Ar/fcaT) % - A r ). (44) 

As expected, the single 5-peak splits into two peaks but there is no broadening. The 
peak separation is |A + — A_|, see (|3~%j). and the peak heights are generally different. To 
lowest order in X/B, the two peaks are located symmetrically around uj = B. Notice 
that for B perpendicular to the effective SO vector A, the splitting is zero. 

It should be stressed that these results hold both for the non-interacting and the 
interacting case. However, for the interacting case realized in SWNTs the double peak 
structure is only possible if spin-charge separation is present. Otherwise the charge 
sector will mix in, leading to broad bands with threshold behaviours ^3]. Closer 
inspection shows that inclusion of the subleading term (J57)l preserves the splitting into 
two peaks, but the peaks now acquire a finite width ~ |A'|. Similarly, the effects 
of backscattering are expected to be small. It is important to point out that the 
double peak spectrum would also show up in the non-interacting case, and hence its 
experimental observation does not represent a true proof for spin-charge separation. 
However, it would certainly give strong evidence, as in the interacting situation the 
double peak survives only when spin-charge separation is realized. 

In practice, to get measurable intensities, one may have to work with an ensemble of 
SWNTs. The proposed experiment may be possible using electric-field-aligned SWNTs, 
or by employing arrays of identical SWNTS. In more conventional samples with many 
SWNTs, however, the SO vector A can take a random direction. For results on the 
averaged spectrum, see [TU] . 

5. ESR spectra for armchair SWNTs: Two-leg Hubbard ladders 

In this section we discuss an alternative description of the ESR spectrum via a lattice 
fermion description. We shall numerically compute ESR spectra for a Hubbard chain 
and a two-leg Hubbard ladder jHj. These models ignore the long-range part of the 
Coulomb interactions which makes them less accurate for the description of charge 
transport. Since the long-range tails of the interaction do not affect the spin dynamics, 
however, this does not create a problem for the spin problem at hand. 

Modelling the SWNT by a two-leg Hubbard ladder is attractive for the following 
reasons: 

• Backscattering interactions as well as commensurability effects (Umklapp 
scattering) are included. 

• The two flavour degrees of freedom are taken care of. 
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• Band curvature effects are taken into account. 

• The marginal SO term (|37j) as well as all irrelevant contributions can be kept. 

• For the Hubbard chain, it is known that any magnetic field spoils spin-charge 
separation |18j . The importance of this effect for SWNTs can be assessed. 

• Accurate numerical techniques are available for such models. 

To gain a broader perspective and to also make contact to other ID conductors, before 
dealing with the more complex two-leg ladder system, we first consider a ID Hubbard 
chain jU] with uniform SO coupling, 

5.1. Hubbard chain 

— * 

The ID Hubbard model with a uniform SO coupling vector A is defined by 

H = £ [4 ("* + c n+ i + H.c] n n] n nl - B £ S z n . (45) 

n n n 

For A = 0, equation (|3H|) leads to a LL phase at low energies It is worth noting that 
then the exact solution of ()45|) shows that with a magnetic field, the charge and spin 
sectors get mixed [18] . To make progress, we gauge away the SO term by the unitary 
transformation 

( c n (n odd) , 

[ exp(ia A • a) c n (n even), 



d r , 



where a = arctan(A/t), A = A/A with A = |A|. The hopping amplitude then changes as 
t — > t = t + A sin a, and one can easily see that n n ^n n ^ = n n ^n n i, where h nu = d^d^. 
Then the Hamiltonian (|43|) can be rewritten as 

H = -t^2(dld n+1 + H.c.) + UJ2 n n1 h ni -BY<S z n , 



Sn 



and the spin operator as 

J \d\pdn (n odd) , 

\ ldl(e taX -°ae- iaX -°)d n (neven). 

The transformed picture is easier to use in practical calculations. We have calculated 
the ESR spectrum for ()45|) using exact diagonalization for small lattices. 

The main features emerging from exact diagonalization can be seen in figure [0 
For weak interactions, see figure ^a), we recover the two-peak structure predicted by 
field theory, see (J44|) . with the peaks symmetrically arranged around uj = B but with 
different peak heights. In addition, no splitting is seen for B perpendicular to the SO 
vector A. The interaction dependence of the ESR spectrum is shown in figure ^b). 
Strong interactions suppress one of the peaks and enhance the other, and also give rise 
to more structure in the spectrum. This behaviour can be traced back to the effects of 
band curvature. Note that for increasing U, the distance between both peaks decreases. 
This feature can be qualitatively explained by a decrease of the spin velocity v s with U, 
i.e. for a flatter band, the spin-orbit splitting is less effective. 
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Figure 1. Exact diagonalization results for the T = ESR spectra of a Hubbard 
chain with A = 0.1, B/t = 0.2, and artificial broadening 0.05i. In (a) ESR spectra for 
U/t = 1 are shown for different angles 7 between A and B. In (b) the [/-dependence 
is shown for 7 = 0. All results are for 12 lattice sites and 4 electrons. 



5.2. Two-leg Hubbard ladder: Armchair SWNTs 

We are now ready to generalize the above treatment to armchair SWNTs. Mapping 
the interacting honeycomb lattice Hamiltonian for armchair SWNTs onto the two-leg 
Hubbard model [Hj, one finds 

#0 = J2(- tc li c (n+i)i ~ ^ c li c «2 + H.c.) + n n ,i^n nM . (46) 

ni n,i 

While under the mapping, t± = t up to 1/R corrections pi] , we shall allow for t± 7^ t to 
have better numerical accuracy. Here we also want to include the spin-orbit interaction 
in this mapping. We use the SO vectors (J3TJ), and for simplicity set u n= \ = 0, 
since this SO vector does not appear in the leading low-energy SO vector A in for 
armchair SWNTs. This simplification also allows for numerically feasible calculations. 
Writing w n = (— l) n W + w, where w = («2 — uz)/2 and W = («2 + U3)/2, the spin-orbit 
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Figure 2. Sketch of the two- leg Hubbard ladder including spin-orbit couplings of the 
type considered here. The arrows indicate the phase for spins up. Spins down have 
the opposite phase. 



term can then be written as 



i(cUw. 



° C n+l,l + Cn/n+l • OX n+12 ) + H.C., 



(47) 



consisting of a uniform and an alternating contribution, 



C n,l/2^n ■ <?C n +l,l/2 ~ C n,l/2 



■ ° C n+l,l/2 ± (- 1 ) n cI,l/2^ • ^Cn+1,1/2- 



The resulting model is studied in the remainder of this section for the special case that 
both w and W are parallel to the uniform magnetic field. Note that w corresponds to the 
leading SO vector A in the low-energy theory, while W is related to the subleading vector 
A', see equation ()36|) . Moreover, away from half-filling, following standard reasoning, 
one may expect that the alternating terms average out. Nevertheless, they will be kept 
below, but we indeed confirm that in large systems they lead only to small effects. When 
the above gauge transformation is applied again, we arrive at a two-leg Hubbard ladder 
(jUIJ), but with the on-chain hoppings t carrying both a uniform phase exp(z0 u ) and an 
alternating phase exp(z0 a ) determined by w and W, respectively, see figure 121 

Let us first discuss the noninteracting case. For U = 0, the solution of equations (jlEJ) 
and (|47|) is straightforward. The eigenenergies are given by (a is the lattice constant) 



E l± {k x ) 
E 2± {k x ) 



(2t cos(k x a + <70 u )) 2 ± 4tt± cos(k x a + cr0 u ) cos( 



t' 



'(2t cos(k x a + 4> u )) 2 ± 4ttj_ cos(k x a + <p u ) cos(^ a ) + t 2 ± 

where a = ±1 for up or down electrons respectively. In figure El we show the non- 
interacting T = ESR spectrum in the presence of both the uniform and the alternating 
phase. The main excitations contributing to each peak are sketched in figure Efb). The 
uniform phase case (see dotted line) produces a splitting of the principal Zeeman peak. 
Due to the different parity of the k = and k = it bands, for the pure uniform case 
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there is no ESR interband transition. When the alternating phase is turned on, the 
folding of the Brillouin zone produces peaks at higher energies. As the parity symmetry 
is broken, interband transitions are now possible. The opening of the small gap due 
to the alternating phase is manifested in the splitting seen in peaks A, B and C. The 
broad structure (F and D) between the highest peaks correspond to transitions between 
branches with opposite curvature. 

Next we discuss results for U > using the density matrix renormalization group 
(DMRG) technique [THll^D!- We have used periodic boundary conditions keeping 256 
states and the finite-system algorithm. The ESR spectrum (J23j) at T = has been 
calculated using the dynamical DMRG technique [2U 122] for a two-leg ladder with 
uniform phase 4> u only. Additional small alternating phases a did not change the 
spectrum significantly. In figure H] we present these results for a quarter filled chain 
and t± = 1.2 and t — 1. These parameters are used in order to clearly see the effect 
of the SO coupling, since for this filling the excitations occur near the linear term of 



Interacting electrons in nanotubes 



18 







4 rungs 






U=0 






U=1 


J 


L i 


— . — 


1 





1 2 



20- 



15- 



10- 





12 rungs 






U=0 






U=1 


J 




V 



20. 



15- 



10- 







I 1 ' 

8 rungs 

U=0 
U=1 




— J - 




1 





1 ; 






i i ■ 

16 rungs 
U=0 
U=1 


J 


V J 





(co-B)/t 



(co-B)/t 



Figure 4. ESR spectra for periodic quarter-filled ladders with a uniform phase 
(j) u = 0.1, t = 1 and t± = 1.2. 



the cosine bands for system sizes that are multiple of 4. The effect of the correlations 
mainly consists of a small shift of the peaks, but the double peak spectrum is preserved. 
Therefore the numerical results lend support to the basic prediction of the analytical 
low-energy theory, and show that the expected double-peak ESR spectrum is stable 
with respect to the above-mentioned perturbations. 

6. Conclusions 

We have reviewed the analysis of the ESR spectrum produced by the spin-orbit coupling 
in SWNTs. The effective field theory analysis shows that at low energy the SO 
interaction only acts in the spin sector and the single Zeeman peak, characteristic of a 
system with SU (2) spin symmetry, splits into two peaks with no broadening. This result 
relies in an essential way on the property of spin-charge separation characteristic of the 
Luttinger liquid state realized in the SWNT. Thus the observation of such a splitting 
would point to this elusive feature. 

The field theory analysis has been complemented by the exact numerical calculation 
of the ESR spectrum for Hubbard models on a single chain and on a two-leg ladder. 
While these models are not able to describe reliably the charge transport properties 
of SWNTs due to neglect of the long-range tails of the Coulomb interaction, they 
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are appropriate descriptions of this spin problem. The numerical analysis confirms 
the field theory predictions, including the validity of the approximations involved in 
their derivation. In addition, it reveals additional structure in the spectrum due to 
band curvature and higher energy processes, which are not captured by the field theory 
approach. 
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